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Abstract 

Noise-enhanced chaos in a doped, weakly coupled GaAs/Alo.45Gao,55As superlattice has been 
observed at room temperature in experiments as well as in the results of the simulation of nonlinear 
transport based on a discrete tunneling model. When external noise is added, both the measured 
and simulated current-versus-time traces contain irregularly spaced spikes for particular applied 
voltages, which separate a regime of periodic current oscillations from a region of no current 
oscillations at all. In the voltage region without current oscillations, the electric field profile consist 
of a low-field domain near the emitter contact separated by a domain wall consisting of a charge 
accumulation layer from a high-field regime closer to the collector contact. With increasing noise 
amplitude, spontaneous chaotic current oscillations appear over a wider bias voltage range. For 
these bias voltages, the domain boundary between the two electric-field domains becomes unstable, 
and very small current or voltage fluctuations can trigger the domain boundary to move toward the 
collector and induce chaotic current spikes. The experimentally observed features are qualitatively 
very well reproduced by the simulations. Increased noise can consequently enhance chaotic current 


oscillations in semiconductor superlattices. 


I. INTRODUCTION 


With the continuous advancement of nonlinear science and the theory of random dynami- 
cal system, we have become aware of the fact that noise should not always be considered as a 
disturbing factor. Actually, a small amount of noise may enhance the dynamics of a system 
so that it becomes better defined and controllable. Constructive effects of noise in nonlinear 
systems have been investigated extensively in the context of stochastic resonances and co- 
herence resonances. By stochastic resonance, noise can optimize the response of a system to 
an external signal and induce stochastic phase synchronization to the external force [1, 2]. 
For a coherence resonance, pure noise without an external signal can generate the most 


coherent motion in the system, as has been mainly observed in excitable systems [8, 4]. 


A doped, weakly coupled semiconductor superlattice (SSL) represents an almost ideal 
spatially one-dimensional nonlinear dynamical system with a large number of degrees of 
freedom, the nonlinearity of which is due to sequential resonant tunneling between adjacent 
quantum wells. Fluctuations of the layer thicknesses, electron density, energy levels, and 
inter-well coupling transform a weakly coupled SSL into a complex nonlinear system, in 
which the electron transport is strongly dissipative. A great richness of nonlinear transport 
behavior has been observed in weakly coupled SSLs, including the formation of stationary 
electric-field domains, periodic as well as quasi-period current self-oscillations, and even 
driven as well as undriven chaos [5, 6]. The oscillatory behavior is attributed to the localized, 
oscillatory motion of the domain boundary, which separates the high from the low electric- 
field domain. Only very recently, spontaneous chaotic [7] and quasi-periodic [8] current 
self-oscillations were observed at room temperature in GaAs/(Al,Ga)As SLs using an Al 


content of 4596, which results in the largest direct barrier for this materials system. 


Both theoretical and experimental results have proven that noise can affect the charge 
transport in weakly coupled GaAs/AIAs SLs [9, 10]. A noise-enhanced coherence resonance 
has been observed in a weakly coupled GaAs/A1As SL at 77 K [11]. The main reason for the 
existence of such a coherence resonance in such a SL is related to the interaction between the 
noise and the two oscillation modes existing in weakly coupled SSL oscillators such as the 
well-to-well hopping mode and the dipole-motion mode. On this basis, numerical simulations 


have shown that noise enhances spontaneous chaos at room temperature [12, 13]. 


In this paper, we report clear evidence for noise-enhanced spontaneous chaos in a doped, 


weakly coupled GaAs/Alo.45Gao.55As SL at room temperature. We find that, with increasing 
noise amplitude, the spontaneous chaotic oscillations appear over a wider range of voltages. 
The experimentally observed results are qualitatively very well validated by numerical sim- 
ulations based on a discrete resonant tunneling model that captures the main features of 
vertical transport in doped, weakly coupled SSLs [12, 13]. The paper is organized as fol- 
lows. In the next section, we describe the sample structure and measurement techniques 
followed by the experimental results of the current-voltage characteristics, the current self- 
oscillations, and attractors extracted from the experimental results. In Sec. III, we present 
the discrete resonant-tunneling model and the results of the numerical simulations. Finally, 


we summarize the obtained results and conclude in Sec. IV. 


Il. EXPERIMENTAL RESULTS 


A. Sample structure and measurement techniques 


The sample consists of a doped, weakly coupled GaAs/(Al,Ga)As SL with 50 periods, 
each period consisting of a GaAs quantum well and an Alo45Gao.55As barrier. The central 
3 nm of each 7-nm-thick GaAs well are doped with Si at a density of 2x10!" cm ?. The 
thickness of the Alo.45Ga9.55AÀs barriers is 4 nm, resulting in a rather weak coupling between 
adjacent quantum wells. The SSL is sandwiched between two highly doped GaAs contact 
layers forming an n*-n-n* diode. For more details of the sample structure, see Ref. 14. 
After plasma etching and providing Ohmic contacts of AuGe/Ni/Au with 35/10/500 nm us- 
ing electron beam evaporation followed by rapid thermal annealing at 420 °C, square mesas 
with a side length of 30 um are investigated. Single devices are wire bonded for electrical 
measurements. All experimental measurements are performed at room temperature. The 
DC bias voltage Vpc and the noise with amplitudes Vise are applied using the function 
generator Agilent 33220A. The current-time traces and the frequency spectra are recorded 
with a 6 GHz oscilloscope LeCroy Wavepro 760Zi-A. For the current-time traces, the os- 
cilloscope records voltage pulses, which are then converted in current values using the /-V 
characteristics. A schematic of the measurement circuit is shown in Fig. 1, which is the 


same as the one in Ref. 11. 
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FIG. 1. Schematic of the experimental setup. The DC bias voltage Vpc and the noise with 
amplitudes Vnoise are supplied by a function generator. The current-time traces and the frequency 


spectra are recorded with a 6 GHz oscilloscope over a 50 €) resistor. 
B. Current-voltage characteristics 


A typical current-voltage characteristic of the the doped, weakly coupled SSL device is 
shown in Fig. 2. In the voltage range between 850 and 860 mV, the DC current rapidly 
increases from 1.24 to 1.55 mA, followed by a current plateau between 860 and 930 mV. 
In this voltage range, electric-field domains are formed with a charge accumulation layer 


forming the domain boundary. As the voltage increases further to 940 mV, the current 
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FIG. 2. Measured current Ipc versus voltage Vpc of the doped, weakly coupled 
GaAs/Alo.45Gao.55 As SL at room temperature. At Vpc = 840 and 940 mV, the frequency spectra 


are dominated by broad-band noise. 


decreases to 1.50 mA and then continues down to 1.31 mA. 


C. Current self-oscillations and frequency spectra 


The frequency spectra displayed in Fig. 3(a) recorded between Vpc = 800 and 960 mV 
without any external noise do not show any current self-oscillations for Vpc values between 
840 and 940 mV. For Vpc values outside this region, i. e., smaller than 840 or larger than 
940 mV, the frequency spectra exhibit periodic oscillations with higher-order harmonics, and 
the fundamental frequency depends on Vpc. At Vpo = 840 and 940 mV, where the transi- 
tions between current oscillations and no current oscillations occur, the frequency spectra 
are dominated by an abrupt decrease to zero frequency. These transitions are accompanied 
by an abrupt change of the current as shown in Fig. 2. This observation may indicate 
an infinite period collision of the oscillatory attractor with a homoclinic orbit. Numerical 


simulations confirm this scenario for the SSL model. 
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FIG. 3. Measured frequency spectra of the current versus Vpc recorded between 800 and 960 mV 


(a) without noise, (b) with Vnoise = 20 mV, and (c) with Voi = 80 mV. 


Figures 3(b) and 3(c) show the frequency spectra for applying additional external noise 
with amplitudes of Vacise = 20 and 80 mV, respectively. For V444, = 20 mV in Fig. 3(b), 
the chaotic oscillations in the transition regions appear over a larger voltage interval than 
without any external noise |cf. Fig. 3(a)]. The voltage region with no current oscillations 
becomes narrower. Increasing the noise amplitude to 80 mV results in a further shrinkage of 
the voltage region with no current oscillations, and the voltage range of seemingly chaotic 
oscillations extends to 855 and 929 mV as shown in Fig. 3(c). At the same time, the frequency 
spectra become broader with increasing noise amplitude. These experimental data clearly 
demonstrate that noise can significantly enhance the chaos generated in doped, weakly 
coupled SSLs. This was shown theoretically for the case of internal (shot and thermal) noise 
by the results of the numerical simulations in Ref. 12. The present results show that chaos 


enhancement can also be achieved more generally by adding appropriate external noise. 


Current traces as a function of time have been recorded for two noise amplitudes at 


Voc = 939 mV as shown in Fig. 4. Note that the current traces are recorded only with AC 
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FIG. 4. Measured current oscillations J (t) recorded for different noise amplitudes (a) Vaoise = 20 mV 


and (b) Vaois = 80 mV. Vpc was fixed at 939 mV. 


coupling so that the information about the DC current level is lost. For Vnoise = 20 mV, 
the current trace contains several irregularly spaced spikes as shown in Fig. 4(a). When 
the noise amplitude is increased to 80 mV, the current trace exhibits more spikes with 
somewhat smaller amplitudes as shown in Fig. 4(b). Similar current traces are observed for 
Voc = 844 mV (not shown). There are typically two oscillation modes present in doped, 
weakly coupled SSLs, the dipole motion mode and the well-to-well hopping mode. When 
the noise level is low, the dipoles are formed at the emitter moving toward the collector with 
a large interval of noise spikes. As the noise amplitude is increased, a continuous motion of 
the well-to-well hopping of the domain boundary occurs, introducing fast current oscillations 
in the SSL. 

The power spectrum shown in Fig. 5(a) recorded at Voc = 944 mV clearly demonstrates 
that, even without any external noise, the doped, weakly coupled SSL can already exhibit a 
broad-band spectrum, indicating the possibility of the presence of chaotic oscillations even 
without external noise. When noise is added as in Fig. 5(b), the power spectrum becomes 


even broader. At the transition regions between periodic oscillations and no oscillations, 


(a) without noise 


(b) VV -20mV 


noise 


Power (dBm) 


0.0 0.5 1.0 1.5 2.0 
Frequency (GHz) 


FIG. 5. Measured frequency spectra recorded (a) without any external noise and (b) with Vhoise = 


20 mV. Vpc was fixed at 944 mV. 


the system appears with increasing external noise amplitude to exhibit a chaotic attractor, 
which leads to spontaneous chaos at room temperature. As predicted by the results of the 
numerical simulations reported in Ref. 12, internal noise enhances spontaneous chaos. The 
same effect of chaos enhancement is observed when we add external noise with increasing 
amplitude as will be further discussed below. At the same time, the frequency spectra in the 
periodic region are dominated by an oscillation with a single frequency, but with increasing 
bandwidth as the noise amplitude is increased. Noise can significantly enhance the chaos 


generated in doped, weakly coupled SSLs. 


D. Attractors obtained from experimental results 


Figures 6(a), 6(b), 6(c), and 6(d) show reconstructed Poincaré maps from the experimen- 
tal current traces versus time for no noise as well as Vacise = 20, 80, and 120 mV, respectively, 
at Voc = 845 mV. In the absence of external noise, the deterministic attractor is station- 
ary, and the Poincaré map should consist of a single point. In Fig. 6(a), the unavoidable 
internal noise produces many points concentrated in a small region. When external noise 
is added, the Poincaré maps explore a region that becomes larger and fuller as the external 
noise amplitude increases as shown in Figs. 6(b), 6(c), and 6(d). Similar Poincaré maps 
are extracted from experimental current traces versus time at Voc = 940 mV. Figure 6(e) 
shows the values of the capacity dimension Do, the information dimension D, and the cor- 
relation dimension Də as a function of the external noise amplitude, which all increase with 
increasing noise amplitude. Without external noise, the SSL is in a stationary state with 
zero fractal dimension. When the external noise is turned on and increased, time-dependent 
oscillations appear, and the fractal dimension Dy increases abruptly demonstrating evidence 
of chaotic oscillations induced by noise. Figure 6(f) shows the multifractal dimension for 
the different noise amplitudes. For large positive values of q, the fractal dimension D, cor- 
responds to regions of the attractor that are more often visited by the system trajectory 
and does not increase that much with external noise. For q < 0 and |q| large, D, is due to 
the least sampled regions of the attractor, and it increases more strongly when the external 
noise increases. We conclude that increasing the external noise expands the regions in phase 
space that are visited by the system trajectory, an effect that is obvious from the sequence 


shown in Figs. 6(b)-6(d). 
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FIG. 6. Poincaré maps extracted from experimental current traces versus time for (a) no external 
noise as well as with external noise of amplitude (b) Vnoise = 20, (c) 80, and (d) 120 mV. (e) 
Values of the capacity dimension Do (diamonds), the information dimension D, (squares), and 
the correlation dimension D» (circles) as a function of Vnoise- (f) Multifractal spectra. Vaoie = 20 


(solid line), 80 (dashed line) and 120 (dash-dotted line). Vpc was fixed at 845 mV. 
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IIl. MODEL AND SIMULATED RESULTS 
A. Discrete resonant tunneling model 


Numerical results are obtained using a discrete resonant tunneling model that captures 
the main features of doped, weakly coupled SSL [5, 6]. The model consists of dynamical 
equations describing the evolution of the variables such as the electric field — F; as well as 
the two-dimensional (2D) electron density n; at well i, the tunneling current density J; ,;i1 


from well i to i + 1, and the SL total current density J(t): 


dF; 
CHE + Jisi+1 = J(t), (1) 
en; = 
Jiu = l v(E) mni Jig nep), (2) 
* k T _ eFjl ah? n, 
d ea Putin = ap v) In h +e Mt (« maT ES JI , (3) 
€ 
ni = Np- ZF EP (4) 
N 
Voc + n(t) 
Yo = POE ult) = malt) + lt) (5) 
i=1 
n 
Jos1 = Coto,  JNSN41 = oon Fn. (6) 
D 


Here i = 1,..., V, where N denotes the number of periods of the SSL. The forward electron 
velocity v\)(F;) is a function with peaks corresponding to the discrete energy levels in 
every well calculated using a Kronig-Penney model for the investigated SSL configuration 
as summarized in Tab. I. A more detailed description can be found in Refs. 6 and 12. The 
overall voltage drop between the ends of the SSL Xo F;/l is expected to be equal to the 
DC voltage bias Vpc and the noise voltage drop n(t) provided by the function generator in 
Fig. 1. The numerical values of the parameters correspond to the experimental configuration 
described in Sec. II with an equivalent 2D doping density due to the doping of the central 
part of the quantum well of Np = 6x10!° cm~?. The effective mass of the electrons in the 
GaAs/Al,Ga,_,As SL is m* = (0.063 --0.083x)m,, where m, denotes the free-electron mass. 
The well and barrier widths are lẹ = 7 and lj = 4 nm, respectively, as for the experimentally 
investigated SSL so that the period of the SSL is | = lj + lu = 11 nm. The transversal area 
of the SSL is assumed to be A = s? with s = 30 pm. Finally, the relative permittivity of 
the SSL is € = l/[lu/€w + ly/e,] with e, = 12.9eg and ep = 10.9es referring to the relative 
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TABLE I. Values of the potential barrier V54,, and the first three energy levels E41, E2, and Eg for 
the doped, weakly coupled GaAs/Alo.45Gag.55As SL. 


Voarr (meV) | E, (meV) ə (meV) £3 (meV) 


388 45 173 346 


permittivity of the well and barrier material, respectively. eo, —e < 0, kg, T, and oy denote 
the vacuum permittivity, the electron charge, Boltzmann’s constant, the lattice temperature, 
and the contact conductivity, respectively. The contact conductivity oo is derived from a 
linear approximation of the emitter current density Jo_,;, which depends on the structure of 
the contact. We use oo = 0.763 A/(V m), which qualitatively reproduces the experimental 
current-voltage characteristic shown in Fig. 2 for a barrier height Via, = 388 meV. The 
voltage noise 7(t) has two components: (i) 7),(t), which is related to the intrinsic noise 
of the source, and (ii) the external noise 7,(t). The thermal noise a(t) is simulated by 
picking random numbers every 5x10^!! s from a zero mean distribution with a standard 
deviation of 1.6x10^? V. The external noise 7.(t) is simulated by picking random numbers 


every 5x10? s from a zero mean distribution with a tunable standard deviation. 


B. Results of numerical simulations 


Figure 7 depicts the simulated current-voltage characteristic for the deterministic equa- 
tions of our model |n(t) = 0)]. The mean current (dashed line) should be compared with 
the experimental one in Fig. 2. The average current for voltages at which there are current 
oscillations in Fig. 7 is qualitatively similar to the first part of Fig. 2. There is a shift 
between the theoretical and experimental curves due to the built-in voltage and a voltage 
drop at an external resistance that are absent in our model equations. In the experiment, 
there is a voltage region of stationary current between two regions of current oscillations. 
This window of stationary currents is absent in the simulations of our model. However, the 
end of the first voltage region with current oscillations in Fig. 2 is similar to the end of the 
voltage region with current oscillations in Fig. 7. 

Figures 8(a), 8(b), 8(c), and 8(d) show the results of the numerical simulations for the cur- 
rent oscillations for different noise amplitudes (a) Vnoise = 1 (only thermal noise is present), 


(b) 8, (c) 33, and (d) 46 mV, respectively, at Vpc = 380 mV, which lies in the oscillatory 
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FIG. 7. Results of the numerical simulations of the current Ipc versus voltage Vpc for o = 
0.763 A/(V m). The dash-dotted, the dashed, and the solid line correspond to the maximum, 
mean, and minimum current, respectively. In the voltage region between 0.26 and 0.4 V, current 
oscillations occur. The simulated curves are shifted to lower voltages compared with the experi- 
mental curve due to built-in fields and a voltage drop at an external resistance. Nevertheless, the 


oscillatory region is qualitatively similar to the one found in the experiments as shown in Fig. 2. 


region of the current-voltage curve presented in Fig. 7. The internal and external noise are 
numerically generated as explained above. With increasing Vhoise, the oscillations become 
less periodic, contain more spikes in the same time interval, and the amplitude becomes 
more and more irregular. The corresponding frequency spectra as determined by a numer- 
ical Fourier transform are displayed in Fig. 8(e), 8(f), 8(g), and 8(h), respectively. The 
number of peaks in the frequency spectra decreases with increasing noise amplitude, while 
the width of the frequency spectra increases. We have repeated the simulations for a voltage 
of 460 mV, for which the deterministic system is in a stationary state. The behavior with 
increasing noise amplitude is very similar to the one shown in Figs. 8(b)-8(d) with some 
differences in the details of the current oscillations and frequency spectra. In comparison 
to Fig. 8(a), it looks of course different, because there are no current oscillations in this 
case. Overall, the experimental trend with increasing noise amplitude is reproduced, i. e., 
the larger the controlled noise amplitude, the more random the current oscillations become. 
Note that simulations of the field distribution (not shown) indicate that the current spikes 
in the simulated current traces in Figs. 8(a)-8(d) are caused by the generation of a small 


charge dipole, which is generated at the emitter contact and moves toward the collector 
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FIG. 8. Results of the numerical simulations for the current oscillations I(t) for different noise 
amplitude (a) Vnoise = 1, (b) 8, (c) 33, and (d) 46 mV at Vpc = 380 mV in the oscillatory region 
of the current-voltage characteristics presented in Fig. 7. The corresponding frequency spectra as 


determined by a numerical Fourier transform are shown in (e), (f), (g), and (h), respectively. 


contact. When the front part of the dipole reaches the collector contact, it disappears, while 


the back part of the dipole rebuilds the domain wall, which existed before the current spike. 


We have also constructed phase portraits from the calculated electric field in well number 
47 versus well number 3. We can build alternative Poincaré maps based on the intersection 
of the trajectories in these phase portraits with an appropriate segment of a straight line. 


Period doubling scenarios or quasi-periodic attractors could be visualized using such Poincaré 
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maps. Thus, these phase portraits give information on the chaotic state of the system. 
Figures 9(a), 9(b), 9(c), and 9(d) show the phase portraits for different noise amplitude 
(a) Vaoise = 1, (b) 8, (c) 33, and (d) 46 mV, respectively, at Vpo = 380 mV, which lies 
in the oscillatory region of the current-voltage curve presented in Fig. 7 and corresponds 
to the simulated data in Fig. 8. When only thermal noise is present |cf. Fig. 9(a)], the 
oscillations are periodic, even though they are somewhat distorted by noise. However, when 
controlled noise is added as shown in Figs. 9(b)-9(d), the oscillations loose their periodicity 
and become more and more distorted with increasing noise amplitude, i. e. the complexity 
of the phase portraits strongly increases, indicating increasing chaotic behavior. The phase 
portraits for a voltage of 460 mV (not shown), for which the deterministic system is in 
a stationary state, look qualitatively very similar to the ones shown in Fig. 9, except for 
the first one with only thermal noise, because there are no current oscillations in this case. 


Note that the thermal noise does not distort the periodicity, since it only introduces some 
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FIG. 9. Results of the numerical simulations for the phase portraits, where the electric field in well 
number 47 is presented versus the electric field in well number 3. The noise amplitude increases 
from (a) Vnoise = 1, (b) 8, (c) 33, to (d) 46 mV at Vpc = 380 mV in the oscillatory region of the 
current-voltage characteristics presented in Fig. 7. Loops correspond to oscillations depicted in 


Fig. 8. 
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small fluctuations. However, by introducing a longer time-correlated (controlled) noise, clear 


chaotic behavior is obtained. 


IV. SUMMARY AND CONCLUSIONS 


We have observed that the external noise can induce spontaneous chaotic oscillation in 
a narrow voltage interval in a doped, weakly coupled GaAs/AIAs SL at room temperature. 
Results of numerical simulations of nonlinear transport based on a discrete tunneling model 
qualitatively confirm the experimentally observed features. While in a noise-free SSL static 
domain boundaries are formed, the domain boundary moves into the collector, new do- 
main walls are formed near the emitter, and chaotic current oscillations are induced, when 
noise is added. Therefore, with increasing noise amplitude, chaotic current oscillations can 
be enhanced in semiconductor superlattices. This approach is consequently a very robust 
method of producing chaotic current oscillations in doped, weakly coupled semiconductor 


superlattices. 
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